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TECHNOLOGY OPERATIONS 


The Aerospace Corporation functions as an "architect-engineer" for national security 
programs, specializing in advanced military space systems. The Corporation's Technology 
Operations supports the effective and timely development and operation of national security 
systems through scientific research and the application of advanced technology. Vital to the 
success of the Corporation is the technical staffs wide-ranging expertise and its ability to stay 
abreast of new technological developments and program support issues associated with rapidly 
evolving space systems. Contributing capabilities are provided by these individual Technology 
Centers: 

Electronics Technology Center: Microelectronics, VLSI reliability, failure 
analysis, solid-state device physics, compound semiconductors, radiation effects, 
infrared and CCD detector devices, Micro-Electro-Mechanical Systems (MEMS), 
and data storage and display technologies, lasers and electro-optics, solid state laser 
design, micro-optics, optical communications, and fiber optic sensors; atomic 
frequency standards, applied laser spectroscopy, laser chemistry, atmospheric 
propagation and beam control, LIDAR/LADAR remote sensing, solar cell and array 
testing and evaluation, battery electrochemistry, battery testing and evaluation. 

Mechanics and Materials Technology Center: Evaluation and characterization of 
new materials: metals, alloys, ceramics, polymers and their composites, and new 
forms of carbon; development and analysis of thin films and deposition techniques; 
nondestructive evaluation, component failure analysis and reliability; fracture 
mechanics and stress corrosion; development and evaluation of hardened 
components; analysis and evaluation of materials at cryogenic and elevated 
temperatures; launch vehicle and reentry fluid mechanics, heat transfer and flight 
dynamics; chemical and electric propulsion; spacecraft structural mechanics, 
spacecraft survivability and vulnerability assessment; contamination, thermal and 
structural control; high temperature thermomechanics, gas kinetics and radiation; 
lubrication and surface phenomena. 

Space and Environment Technology Center: Magnetosphenc, auroral and cosmic 
ray physics, wave-particle interactions, magnetosphenc plasma waves; atmosphenc 
and ionospheric physics, density and composition of the upper atmosphere, remote 
sensing using atmospheric radiation; solar physics, infrared astronomy, infrared 
signature analysis; effects of solar activity, magnetic storms and nuclear explosions 
on the earth's atmosphere, ionosphere and magnetosphere; effects of electromagnetic 
and particulate radiations on space systems, space instrumentation; propellant 
chemistry, chemical dynamics, environmental chemistry, trace detection; 
atmosphenc chemical reactions, atmosphenc optics, light scattenng, state-specific 
chemical reactions and radiative signatures of missile plumes, and sensor out-of- 
field-of-view rejection. 
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Stormtime Ring Current and Radiation Belt Ion Transport: 
Simulations and Interpretations 


Margaret W Chen, Michael Schulz 1 , Larry R Lyons, and David J Gomey 

Space and Environment Technology Center, The Aerospace Corporation, El Segundo, California 


We use a dynamical guiding-center model to investigate the stormtime transport of ring current 
and radiation-belt ions We trace the motion of representative ions’ guiding centers in response to 
model substorm-associated impulses in the convection electric field for a range of ion energies Our 
simple magnetospheric model allows us to compare our numerical results quantitatively with ana- 
lytical descriptions of particle transport, (e g , with the quasilinear theory of radial diffusion) We 
find that 10-145-keV ions gain access to L ~ 3, where they can form the stormtime ring current, 
mainly from outside the (trapping) region in which particles execute closed drift paths Conversely, 
the transport of higher-energy ions ( it 145 keV at L ~ 3) turns out to resemble radial diffusion 
The quasilinear diffusion coefficient calculated for our model storm does not vary smoothly with 
particle energy, since our impulses occur at specific (although randomly determined) times De- 
spite the spectral irregularity, quasilinear theory provides a surprisingly accurate description of 
the transport process for t 145-keV ions, even for the case of an individual storm For 4 different 
realizations of our model storm, the geometric mean discrepancies between diffusion coefficients 
obtained from the simulations and the quasilinear diffusion coefficient amount to fac- 
tors of 2 3, 2 3, 1 5, and 3 0, respectively We have found that these discrepancies between 
and can be reduced slightly by invoking drift-resonance broadening to smooth out the sharp 
minima and maxima in D q L , The mean of the remaining discrepancies between and D£, 

for the 4 different storms then amount to factors of 1 9, 2 1, 1 5, and 2 7, respectively We find 
even better agreement when we reduce the impulse amplitudes systematically m a given model 
storm (e g , reduction of all the impulse amplitudes by half reduces the discrepancy factor by at 
least its square root) and also when we average our results over an ensemble of 20 model storms 
(agreement is within a factor of 1 2 without impulse-amplitude reduction) We use our simula- 
tion results also to map phase-space densities / in accordance with Liouville’s theorem We find 
that the stormtime transport of > 145-keV ions produces little change in / the drift-averaged 
phase-space density on any drift shell of interest However, the stormtime transport produces a 
major enhancement from the pre-storm phase-space density at energies ~ 30-145 keV, which are 
representative of the stormtime ring current 


l Introduction 

This work is an outgrowth of our continuing study 
of energetic charged-particle transport m the magne- 
tosphere. The study began as an effort to understand 
the development of the stormtime ring current but has 
expanded to include the radial diffusion of radiation- 
belt ions Our study entails a guiding-center simula- 
tion of particle motion in the presence of a succession 
of substorm-associated impulsive enhancements in the 
magnetospheric convection electric field We synthe- 
size our model storms by means of a random-number 

x Now at Lockheed Research Laboratory, Palo Alto, California 

Space Plasmas Coupling Between Small 
and Medium Scale Processes 
Geophysical Monograph 86 

Copyright 1995 by the American Geophysical Union 


generator and apply them to a simple magnetospheric 
model m order to make the storm effects realistic but 
mathematically analyzable 

We have found, in agreement with Lyons and 
Williams [1980], that the access of ~ 10-145 keV ions 
(having first adiabatic invariants fj, ~ 1-13 MeV/G) to 
the region (L ~ 3) where they can form the storm- 
time ring current occurs largely as a consequence of 
the enhanced mean- value of the convection electric field 
rather than from its impulsive character Indeed, most 
of the particles in this energy range that reach L ~ 3 
in our model storm turn out to have been transported 
there from outside the (trapping) region in which par- 
ticles execute closed drift paths Conversely, the trans- 
port of higher-energy particles (having, for example, £ 
145 keV at L ~ 3) turns out to resemble radial diffusion 
across closed drift paths [cf Lyons and Schulz, 1989] 
By having formulated the model storm in an eas- 
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lly analyzed way, we are able to compare the radial 
diffusion coefficients obtained from our guiding-center 
simulation with the predictions of quasilinear theory 
[e g , Falthammar, 1965, Cornwall, 1968] and various 
refinements thereof We find that the quasilinear diffu- 
sion coefficient calculated for any of our model storms 
shows a remarkably unsmooth variation with particle 
energy because the impulses occur at specific (although 
randomly determined) times Despite this, quasilinear 
theory provides a surprisingly accurate description of 
the transport process for > 145-keV ions, even for the 
case of an individual storm As expected, the agreement 
becomes even better when we reduce the impulse am- 
plitudes systematically in a given model storm, and also 
when we average our results over an ensemble of model 
storms constructed by the same (random) method 

Of course, a radial diffusion coefficient is not de- 
fined for transport from open to closed dnft trajecto- 
ries, such as we have found to occur for ions having 30 
keV < E £ 145 keV For these energies, however, we 
find a major increase in the drift-averaged phase space 
density / from its pre-storm value upon mapping / in 
accordance with Liouville’s theorem Particle energies 
30-150 keV have been shown by many observational 
studies [e g , Prank, 1967, Smith and Hoffman, 1974, 
Williams and Lyons, 1974, Lyons and Williams, 1976, 
Hamilton et aL, 1988] to be representative of the storm- 
time ring current as a whole In contrast, we find little 
change m drift-averaged phase-space density as a con- 
sequence of stormtime transport for ions having y. > 13 
MeV/G ( E £ 145 keV at R = 3), for which the trans- 
port is diffusive 

2 Field Model 

The magnetic field model that we use m this study 
is obtained by adding a uniform southward field AB 
to the geomagnetic dipole field We invoke this simple 
field configuration because it enables us to make di- 
rect comparisons between the simulated transport and 
previous analytical formulations An advantage of our 
model over a purely dipolar field is the presence of a 
quasi-magnetopause at the boundary between closed 
and open field lines (see Figure 1). The equation of 
a field line m this model is 

[1 + 0 5(r/6) 3 ] -1 (r/I?£;)csc 2 0 = constant = L (1) 

where r is the geocentric distance, 6 is the magnetic 
colatitude, Re is the radius of the Earth, and b = 
1 5 L*Re = 12 82 Re is the radius of the equatorial 
neutral line This value of 6 which is obtained by map- 
ping the last closed field line (denoted L*) to a colati- 
tude of 20° on the Earth, corresponds to |AB| = 14 474 
nT and L* = 8 547 The limit 6 — ► oo (L* — > oo) would 
correspond to a purely dipolar B field In this study, we 



Fig 1 An illustration of the magnetosphenc B-field model used 
in this study The model is symmetric about the sm0 = 0 axis 
and about the equatorial plane, which contains a circular neutral 
line at r = 6 on the magnetic shell L = L* , which approaches an 
asymptotic distance p' from the tail axis at large distances |z| 
from the equatorial plane 

consider only equatonally mirroring particles m which 
the equatorial field intensity Bo is given by 

Bo = ( y E /r Z ) - 14 474 nT, (2) 

where y E = 3 05 x 10 4 nT-Bi- is the geomagnetic dipole 
moment Further details of this field model are given 
by Schulz [1991, pp 98-110] 

We assume that the total electric field E = - Vi*# is 
derivable from the scalar potential 

*— *♦$(£)-*♦ W)-* <*> 

m which the three separate terms correspond to coro- 
tation (Vh — 90 kV), the Volland-Stem [Volland, 1973, 
Stem, 1975] model of quiescent convection (Vo = 50 
kV), and the time-dependent enhancement A V(t) as- 
sociated with the stormtime convection, respectively 
The time-varying term m the potential is assumed to 
vary as L [cf Nishxda, 1966, Bnce, 1967] rather than 
as L 2 because electric disturbances are expected to be 
less well shielded than steady-state convection by the 
inner magnetosphere 

We model the storm-associated enhancement AV(t) 
in the cross-tail potential drop, 

N 

AV * ex p[(** - t )/ r W - **) ( 4 ) 

*=i 
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where 6(t) is the unit step function (= 1 for t > 0, = 0 
for t < 0), as a superposition of almost randomly oc- 
curring impulses that rise sharply and decay exponen- 
tially with a “lifetime” r = 20 mm [cf Cornwall, 1968] 
The impulses represent the constituent substorms of a 
storm The potential drop AV t associated with any in- 
dividual impulse is chosen randomly from a Gaussian 
distribution with a 200-kV mean and a 50-kV standard 
deviation We have chosen such a large mean value of 
AV^ since our intention is to model a major (|D«t| ~ 
200 nT) storm, such as those which Lyons and Williams 
[1980] analyzed Since those storms had a main phase 
lasting ~ 3 hr, we assume that the N start times t, in 
(4) are randomly distributed within a 3-hr time inter- 
val corresponding to the mam phase of a storm How- 
ever, we impose a 10-min “dead time” (after each im- 
pulse onset) during which no subsequent impulse can 
start This constraint imposes a seemingly realistic de- 
lay between consecutive impulse onsets Without such 
a dead-time it would be possible for the next impulse 
to start immediately after the previous one, and this 
could lead to the build-up of unrealistically large cross- 
tail potentials Further details of this model storm are 
given m [ Chen et aL, 19926] 

We have constructed 100 such random storms so that 
on average there are 9 impulses per storm or 3 sub- 
storms/hr We have done this by generating 1800 ran- 
dom numbers between 0 and 100 and disqualifying 
about half of these through the dead-time constraint 
We have randomly chosen four model storms for de- 
tailed case studies Figure 2 shows the variation in 
cross-tail potential for these prototypical storms The 
mean enhancement in cross-tail potential drop for these 
particular storms over the time interval t\ < t < t\ + 3 
hr are (A V{t)) = 180 kV, 178 kV, 154 kV, and 207 
kV, respectively Since we choose to average over the 
period t\ to ti + 3 hr, we may be excluding a signifi- 
cant portion of the last impulse Thus, defined in this 
way, the average cross-tail potential drops are typically 
somewhat less than 200 kV 

3 Particle Dynamics 

Since we simulate the gmdmg-center motion of non- 
relativistic equatonally mirroring particles, we treat the 
first two adiabatic invariants (/i^0 and J = 0, respec- 
tively) as conserved quantities It follows from (2)-(4) 
that the gmdmg-center motion of an equatonally mir- 
roring particle subject to E x B and gradient- 2? drifts 
is described by 


dL 

dt 


L 2 cos <j> 
2 heR-E 2 


v 0 (i) 2 + AK ( o(i) 


(5) 


d<t> _ n „ MMS 

dt ~ “ * qB 0 r* 

where fl is the angular velocity of the Earth, q is the 
charge of the particle, and <f> is the azimuthal coordinate 
(local time) 

We solve the ordinary differential equations (5) and 
(6) simultaneously by using the Buhrsh-Stoer extrapo- 
lation method with variable time-step [e g , Press et aL, 
1986, pp 563-568] for specified initial conditions First 
we obtain steady-state adiabatic drift paths associated 
with a particular value of the first adiabatic invariant 
fi by setting A V(t) = 0 m (5) and (6) Next, to study 
the effects of a time-dependent cross-tail potential, we 
start representative particles at points equally spaced 
m time on a steady-state drift path We do this for 
the purpose of properly calculating a radial diffusion 
coefficient and for obtaining drift-averaged phase space 
densities (see below) (We note that this method dif- 
fers from previous simulation studies [e g , Smith et aL, 
1979, Lee et aL, 1983, Takahasht, 1990] in which parti- 
cles are injected from a mghtside boundary ) We then 
apply (4), which prescribes a storm-associated variation 
A V(t) in the cross-tail potential drop and run the sim- 
ulation to determine the consequent stormtime particle 
transport We can run the simulation either backward 
in time (to determine where any representative particle 
must have been prior to the storm in order to reach the 
desired phase on its “final” drift shell) or forward in 
time (to follow the dispersal of initially co-drifting par- 
ticles among drift shells during and after the storm ) 

4 Simulated Guiding-Center Trajectories 

In this section we present results of simulated storm- 
time transport of singly charged ions having various n 
values and vanishing second invariant (J = 0) For 
illustrative purposes, we present only the results that 
were obtained when we applied the model storm 
((AV(t)) = 180 kV) shown m Figure 2a The dashed 
outer circle on each particle-trajectory plot shows the 
location of the neutral line, a circle of radius R = 12 82, 
which marks the boundary between open and closed 
magnetic field lines in our magnetic-field model (cf 
Figure 1) Figure 3a illustrates steady-state trajecto- 
ries of equatonally mirroring ions for p = 3 MeV/G 
This corresponds to an energy of 33 keV at a geocentric 
radial distance r = 3 Re- For this particular y. value, 
an x-type separatnx marks the boundary between open 
and closed drift trajectones We label closed drift shells 
m terms of the dimensionless third adiabatic invariant 
defined by Roederer [1970, p 1078] as 


1 _ 

$bRe 

_ r i 

r d4> 1 

L ~ 

27 the 

271 j 

m\ 


and 
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TIME, hr TIME, hr 


Fig 2 The cross-tail potential V(t) in our model storm consists of a quiescent value Vo (= 50 kV) and a 
superposition of exponentially decaying impulses (decay time r = 20 min) These impulses represent the constituent 
substorms of a storm and start at times that are distributed randomly over a 3-hr time interval, except that we 
impose (after the start of each impulse) a 10-mm “dead time” after the start of each impulse during which no 
subsequent impulse can start Four realizations of our model storm are shown The average enhancements (AV(t)) 
in cross-tail potential drop over the respective 3-hr main phases are (a) 180 kV, (b) 178 kV, (c) 154 kV, and (d) 
207 kV 


where is the magnetic flux enclosed by that drift 
shell and L(4>) denotes the field-line label at longitude 
<f> on the drift shell In particular, we denote by Li the 
drift shell that separates open and closed drift paths. 

We have examined the access of 10 ns to the quiescent 
drift shell that intersects the dusk meridian at R = 
r/R e — 3 for a range of fx values (3 MeV/G ;$ /x £ 
200 MeV/G) We select this drift shell because it is 
representative of where particles need to be transported 
in order to form the stormtime ring current We find 
the quiescent iomc drift period r 3 to be 12 hr, which is 
longer than the assumed 3-hr main phase Starting with 
12 representative 10 ns (as indicated by the filled circles) 
equally spaced in time on this drift path of interest, we 
have run the simulation backward in time 

The time-reversed trajectories shown m Figure 36 
thus indicate where the particles must have been prior 
to the storm in order to have readied the final drift 
shell (L « 3) of interest We found that 3 of the 12 rep- 
resentative 10 ns would have been transported outward 
from dosed drift paths having smaller L values The 


other 9 representative 10 ns would have been transported 
inward from the night side along open trajectories to 
populate the final closed drift path Eight of these 9 
representative 10 ns would have come from beyond the 
boundary of our model We have compared the time 
(~ 1 — 1 3 hr) required for these 8 representative 10 ns to 
be transported inward from the neutral line to the final 
dosed drift path of interest with the convection time [cf. 
Lyons and Williams, 1980] obtained by integrating (5) 
while keeping cos <j> constant in time Although Lyons 
and Williams [1980] had envisioned direct-convective 
access from dosed drift trajectories, we have found that 
the direct convective access occurs mainly along open 
drift trajectories from the neutral line However, the ac- 
cess times obtained for the simulated trajectories agree 
reasonably well in most cases with estimates based on 
direet-convective access Agreement is especially good 
for 10 ns transported from the mghtside neutral line to 
the vicinity of the final drift path of interest in the 
quadrant centered on midnight Agreement was not as 
good for ions that had drifted to other local times be- 
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p = 10 KeV/G 110 keV 



r, = 2 4 hr (<0 



STEADY STATE 



X(R.) 


Fig 3 Quiet-time equatorial trajectories of singly charged ions having ft = 3, 10 and 30 MeV/G are plotted in the 
upper panels The outer dashed circle represents the neutral line at r = b Ions whose drift paths cross the dusk 
meridian at R = 3 have drift periods T3 as noted The 12 representative ions’ “final” positions on the steady-state 
drift path of interest are denoted by small “filled” circles Corresponding stormtime trajectories computed in our 
time-reversed simulation are shown in the lower panels 


fore reaching the vicinity of the final drift path during 
the storm Reasonably good agreement was typical for 
lower-energy ions (E ~ 30-100 keV) [Chen et al., 19926] 
and thus confirms convective access as their mode of 
transport 

Equatorial steady-state drift paths for ions having 
/x = 10 MeV/G or energies of 110 keV at R = 3 are 
shown m Figure 3c Again, in this case, an x-type sep- 
aratrix marks the boundary between open and closed 
drift paths The quiescent drift period (T 3 = 2 4 hr) 
for these ions on a drift shell that intersects the dusk 
meridian at R = 3 is comparable to the duration of the 
storm’s main phase 

The time-reversed simulated stormtime trajectories 
for ft = 10 MeV/G are shown m Figure 3d At this 
ft value, only half of the representative ions have been 
transported to the final drift shell of mterest by moving 
inward from the night side along open drift trajecto- 
ries The other half have been transported from closed 


drift paths of either smaller or larger L value Their 
stormtime transport begins to resemble radial diffusion 
Thus, the mode of access of intermediate-energy ions to 
the stormtime ring current is transitional between con- 
vective and diffusive 

Finally, we consider ions having ft = 30 MeV/G, 
which would correspond to energies of 335 keV at R = 
3. Steady-state ion drift paths are illustrated in Figure 
3e Here the boundary between open and closed drift 
paths is a closed drift path tangential to the neutral hne 
[cf Brice and Ionanmdts, 1970, Schulz, 1976] Again, 
we consider the access of ions to the drift path that 
intersects the dusk meridian at R = 3 The quiescent 
drift period (73 = 0 73 hr) for 10 ns on such a drift shell 
is less than 1 /4 of the main-phase duration 

The simulated storm transports 10 ns to the drift shell 
of mterest from closed drift paths of either smaller or 
larger L value Thus, there is a spread among the ini- 
tial L values of the representative particles (see Figure 
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3 f] This transport resembles radial diffusion in this 
case, except that this was a time-reversed calculation 
We have also run a simulation forward in time, so as to 
follow the dispersal of ions from a common initial drift 
path Not surprisingly, the resultmg stormtime trans- 
port generates a plot qualitatively similar to Figure 3/ 
although the time-forward simulation applies to an im- 
plicitly different set of particles from those followed m 
the time-reversed simulations Transport of the type il- 
lustrated in Figure 3/ was typical for ions having /i > 13 
MeV/G, which corresponds to energies > 145 keV at 
R= 3 

Figure 3 thus illustrates a range of modes of parti- 
cle access to the stormtime ring current Ions having 
H £ 5 MeV/G (E < 30 keV at R = 3) undergo mainly 
direct convection from open (plasmasheet) drift paths 
to closed drift shells at R ~ 3, while ions having /i £ 13 
MeV/G respond to the same enhancement of the con- 
vection electric field in a manner which resembles ra- 
dial diffusion among closed drift shells [cf Lyons and 
Schulz, 1989] The transition between these two ideal- 
ized modes of access occurs at /i ~ 5-13 MeV/G (E ~ 
55-145 keV at R = 3) for ring-current particles whose 
quiescent drift periods are comparable to the duration 
of the model storm’s main phase 

5 Phase-Space Density Mapping 

We have performed time-reversed ion simulations at 
additional ft values (1 MeV/G £ /t ;$ 100 MeV/G) for 
the purpose of mapping phase space distributions in ac- 
cordance with Liouville’s theorem to drift shells that in- 
tersect the dusk meridian at R = 3 This requires that 
we specify the distribution at the neutral line before 
and during the storm, and the distribution on closed 
trajectories before the storm At the neutral line, we 
maintain an exponential spectrum 

f* = exp(-/i//io) (8) 

for the phase-space density at the neutral line and set 
/io = 5 MeV/G This leads to a reasonable drop off of 
the boundary spectrum at high energies [cf Williams, 
1981] We neglect losses for the distribution along open 
trajectories so that f* specifies the phase space den- 
sity everywhere beyond the boundary between open 
and closed drift trajectories (e g , see upper panels in 
Figure 2), which we label L i(fi) 

We neglect Coulomb drag for simplicity and assume 
that the pre-storm transport of ions on closed trajecto- 
ries is governed by an equation of the form 



where / is the drift-averaged phase-space density at 
fixed /i and J, Dll is the diffusion coefficient for trans- 
port m L, and r q is the ionic lifetime against charge 
exchange The steady-state solution to (9), in which 
radial diffusion balances charge exchange, can be ex- 
pressed m closed form m terms of modified Bessel func- 
tions of fractional order if Dll and T g vary as power 
laws in L [Haerendel, 1968] Thus, we seek to fit Dll 
and Tq accordingly 

We notice, from the plot of H + charge exchange life- 
time profiles (solid curves) reproduced from selected n 
values from Cornwall [1972] in our Figure 4a, that r q 
tends to vary as a power law in L at the smaller L 
values Accordingly, we specify 

T q « L -8 10 3 (/i/4) 4 day -1 (10) 

as a rough approximation corresponding to the dashed 
curves m Figure 4 a This is a fairly good fit to the 
charge exchange lifetimes taken from Cornwall [1972] 
for L ;$ (2 5 /i)w 5 ), which covers most of the range of \i 
and L values of interest However, we hope to improve 
upon our fit of the charge exchange lifetimes in future 
work 



L L 


Fig 4 (a) Profiles of H + charge-exchange lifetimes (solid curves 
taken from Cornwall [1968]) for different values (indicated in 
MeV/G) of the first adiabatic invariant The formula given by 
(10) provides a fairly good power-law fit (L~ 8 , dashed lines) to 
T g for L S (2 5ji)^ 2 ' 5 ) (b) Profiles of Dj\, reciprocal of the 

standard diffusion coefficient given by (11), for selected values of 
H (indicated m MeV/G) Dashed lines represent power-law (L~ 8 ) 
fits specified by reciprocal of (12) 
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Similarly, we estimate a power law fit to the diffusion 
coefficient The standard model [e g , Cornwall , 1972] 
leads to a diffusion coefficient of the form 


Dll 


1 4 x 10- 5 L 10 
/x 2 + L 4 


day \ 


( 11 ) 


where p is in units of MeV/G The solid curves in Figure 
4 b are profiles of D^l f° r selected /i values Because the 
diffusion coefficient varies as L 6 for L 2 3> u and as L 10 
for L 2 <C /x, we have compromised on L* in order to 
obtain a single power law We thus obtam a power-law 
“fit” 

Dll «7x HrV" 1 ^ 8 day -1 (12) 

to the radial diffusion coefficient by requiring that the 
dashed curves m Figure 46 be tangent to the corre- 
sponding solid curves at L = /x 1 / 2 In the future, we 
plan to refine our power-law model for Dll (as well 
as for Tq) However, for the present, our simplistic but 
reasonable power-law fits to the transport coefficients 
allow us to express the pre-storm phase-space distribu- 
tion f((x, L) by means of the equation 


L\ 5/2 /Qx,L) 
Li) /*(a 0 



Fig 5 The pre-storm phase space density f spectrum for ions 
that drift through the dusk meridian at R = 3 is represented by 
the plotted curve The solid and dashed portions of the curve cor- 
respond to closed and open drift trajectories, respectively The 
drift-averaged phase-space density distribution (filled circles) de- 
noted post-mam-phase is obtained by averaging the mapped val- 
ues of / for the 24 representative ions 


-f-5/2(fl)-^-5/2(flo) ~ K- 5/2(0)-f— 5/2(0o) 

/-5/2(0l)^— 5/2(00) ~ ^-5/2(^l K-5/2 (^o). 


, (13) 


where 6 = L(r q DLL )~ ^ and t 9 Dll depends only on 
\i The mner boundary Oq in (13) corresponds to the 
drift shell that grazes the Earth’s atmosphere We ob- 
tam Lq as a very weak function of /x for this purpose by 
evaluating (7) for the drift shell that intersects the dusk 
meridian at R = 1.1 The outer boundary 9\ in (13) 
corresponds to the separatrix L\(n) between closed and 
open drift paths (cf Figure 3, upper panels) 

Using (8) and (13), we plot (see Figure 5) the pre- 
storm phase-space density spectrum /(/x, L ) for the 
drift shell that intersects the dusk meridian at R = 
3 We distinguish between values of /(/x, L) on open 
(dashed curve) and closed (solid curve) drift trajecto- 
ries (/x = 2 7 MeV /G is the smallest first invariant for 
which the trajectories that drift through the dusk me- 
ridian at R = 3 are closed) Our simple model repro- 
duces essential features similar to those found in proton 
phase-space distributions obtained by Wilhams [1981] 
from ISEE 1 data. At the higher /x values for which ra- 
dial diffusion dominates charge exchange, the spectrum 
drops off like our exponential boundary spectrum The 
spectral peak (found at /x ~ 22 MeV/G) had been an- 
ticipated by Spjeldv ik [1977] and occurs mainly because 
the charge-exchange lifetime decreases with decreasing 


/x However, since Li(/x) also varies directly with /x for 
/x > 1 MeV/G, 10 ns havmg small /x values do not have 
to diffuse as far from their boundary between closed 
and open drift paths in order to reach R = 3 (cf upper 
panels in Figure 3) For this reason the solution speci- 
fied by (13) rises again at low /x to join the exponential 
boundary spectrum (dashed curve) which corresponds 
to 10 ns on open drift paths 

We have invoked Liouville’s theorem to map the 
phase-space / for each representative ion from (8) and 
(13) For this portion of the study, we employ time- 
reversed tracings of 24 (cf Figure 3), rather than 12, 
representative 10 ns from end-point phases equally 
spaced in time on a drift shell that intersects the dusk 
meridian at R = 3 By averaging the mapped / values 
for all the 24 representative 10 ns we obtam a good esti- 
mate of the drift-averaged phase-space density / (filled 
circles in Figure 5) attained upon completion of the 
main phase of the model storm Our approach differs 
from that of Ktstier et aL [1989] who made point-to- 
point mappmgs of phase space distributions at vari- 
ous local times using pre-storm spectra obtained from 
AMPTE data. 

We find a major enhancement from the pre-storm 
phase-space for /x ~ 3-13 MeV/G This range corre- 
sponds to energies ~ 30-150 keV, which are known to 
be representative of the stormtime ring current [e g , 
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Lyons and Williams, 1976, Williams, 1981, Ktstler et 
al , 1989] Moreover, these are energies for which our 
simulations have shown that ion transport to L ~ 3 
occurs largely from open trajectories (e g , Figure 3 b) 
In contrast, for p > 13 MeV/G (E > 145 keV at R = 3) 
we find little change m /(m, L) as a consequence of the 
transport associated with a single storm This range 
corresponds to particles whose transport resembles ra- 
dial diffusion (e g , Figure 2 f) 

We find our preliminary results of mapping phase- 
space densities to be particularly satisfying since they 
are consistent with many observed features of ring-cur- 
rent phase-space distributions Thus, we are extending 
this study to other L-shells of interest In addition, we 
are refining our model for the pre-storm and boundary 
phase-space distributions and will report on the results 
in the near future 


6 Diffusion and Quasi-Diffusion Coefficients 


Although the higher-energy ions (E > 145 keV) do 
not seem to contribute much to the stormtime ring cui- 
rent, their diffusive transport is nevertheless interesting 
m the radiation-belt context Thus, we have made de- 
tailed quantitative comparisons of diffusion coefficients 
obtained from our simulations with predictions from 
the quasihnear theory of Falthammar [1968] For this 
purpose, we consider four realizations (shown m Figure 
2) of our random storm model These were randomly 
chosen from the 100 storms that we had originally gen- 
erated 

For each realization of our model storm, we have com- 
puted diffusion coefficients for selected values of p (such 
that 15 MeV/G ^ p £ 200 MeV/G) from the distribu- 
tion of initial L values of the time-reversed trajectories 
(see Figure 3 f) We have done this by constructing the 
quantity 

r\svm — 

U LL = 


X ' L t=l 


12 
L t=i 


■7‘> I 




, (14) 


where L x and L f denote the drift-shell labels of the ini- 
tial and final trajectories (respectively) of the 12 repre- 
sentative ions, and where T (= 3 hr) denotes the dura- 
tion of the mam phase of the model storm The quan- 
tity D” ™ is thus a measure of the variance among the 
initial third adiabatic invariants of particles situated on 
the final drift shell of interest We also computed dif- 
fusion coefficients D"™ from time-forward simulations 
by interchanging the indices i and / in (14) The dif- 
fusion coefficients obtained from simulations run 
forward and backward m time are not very different, 
80% of them being within a factor of 1 6 of each other, 
although they pertain implicitly to different sets of par- 


ticles The geometric mean discrepancies (among the 
12 values of p) between the values of obtained 
from simulations run forward and backward in time 
amounted to factors of 1 3, 16, 1 5, and 1 8 for the 
four model storms in Figure 2 
We compare the diffusion coefficients obtained from 
the simulated trajectories with the resonant^particle 
formulation [Falthammar, 1965] of radial-diffusion the- 
ory in which the diffusion coefficient is of the form 



where E(w/2ir) is the spectral-density of the (quasi- 
uniform) equatorial electric field in the inner magneto- 
sphere and n 3 /27r is the particles’ quiescent drift fre- 
quency When we substitute the spectral-density func- 
tion for our model storm (see Chen et al. [19926] for 
derivation) into (15), we obtain the quasihnear diffusion 
coefficient 


ql t 2 L*R\ * ^ AKAV, cos[fl 3 (t, - *,)] 
LL 167>|(L*)2 l+nfr 2 

(16) 

in which correlations between the impulses lead to cross 
terms (j /i) By neglecting the cross terms, we could 
recover essentially the standard diffusion coefficient of 
Cornwall [1968], but here we retain all the cross terms 
in (16) for comparison with as computed for in- 
dividual storms 

The dashed curves m Figure 6 represent the quasilin- 
ear diffusion coefficients at L « 3 for the corresponding 
four model storms shown in Figure 2 The respective 
quasihnear diffusion coefficients are not very smooth 
functions of p This is because the impulse onsets (i) 
associated with any individual storm modeled by (4) oc- 
cur at specific (although randomly determined) times 
t y , which means that the corresponding spectral density 
E(u)/2tt) is not a very smooth function of frequency 
We plot as data points in Figure 6 the diffusion co- 
efficients D £7* obtained from simulations run forward 
(filled circles) and backward (open circles) in time For 
comparison purposes, we chose soyie of the [i values to 
correspond with the minima m D q ^ L m case (a) Agree- 
ment of the diffusion coefficients obtained from 

the simulations with quasihnear theory js surprisingly 
good despite the strong variability of D g ^ L with fi For 
cases (a)-(d), we find that the geometric means of the 
discrepancies amount to factors of 2 3, 2 3, 1 5, and 3 0, 
respectively Agreement is best for case (c), m which 
there were only 7 substorms during the model 
(see Figure 2c) and consequently less variability l 
with fi We find that quasihnear theory even accounts 
for the n values (e g , p = 75 and 80 MeV/G in Figure 


storm 
i DV, 



hr ' Dll. hr- 
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Fig 6 Plots of diffusion coefficients obtained via (14) from time-reversed (open circles) and time-forward 

(filled circles) simulations, for comparison with the quasilmear diffusion coefficient D V ^ L (dashed curve) given by 
(16) as an implicit function of p for the four realizations of our model storm shown in Figure 2 The diffusion 
coefficients D T £ L (corrected for resonance-broadening) are represented by the solid curves 
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6a) for which the diffusion coefficients computed from 
the simulation are especially small 

However, the diffusion coefficients obtamed from the 
simulated trajectories generally do not show quite as 
much variability with p as quasilinear theory predicts 
This is not surprising, since quasilinear theory postu- 
lates a perfectly sharp resonance at the quiescent drift 
frequency, whereas the simulated transport leads to an 
eventual spread among the quiescent drift frequencies 
of the representative ions for each p A rough estimate 
for the anticipated spread in Q 3 /2ir is 

Aw/2tt « (D LL r/27r 2 ) 1 /2|(aft 3 /aL) M | (17) 

since the mean-square spread m L accumulated during 
the transport is 2DnT [ Chen et ad., 1992a] An es- 
timate for the diffusion coefficient (corrected for reso- 
nance-broadening effects) is thus 


_ L 6 R i 

Du ~ 


rf) 3 +(A»/2) 

Jn 3 -( Au>/ 2 ) V 27I 7 


4/i^Aw 7n 3 -( Au-/2) 


du 


(18) 


Since the frequency bandwidth given by (17) depends 
on Du, we have iterated between (18) and (17) until 
satisfactory convergence to the desired solution (called 
D r ii) is achieved The results for cases (a)-(d) are plot- 
ted as solid curves m Figure 6 We find that inclusion 


of this nonlinear resonance-broadening effect tends to 
reduce the discrepancy between quasilinear theory and 
Dii by a smoothing out the sharp relative maxima 
ana minima with respect to p Corrections were typ- 
ically ~ 10-30% at the relative maxima but were as 
much as 60% at the relative minima (e g , near p=\M 
MeV/G m Figure 6a). The geometric means of the re- 
maining discrepancies between and amount 
to factors of 1 9, 2 1, 1 5, and 2 7 for cases (a)-(d), 
respectively To determine whether the remaining dis- 
crepancies are attnbutable to the neglect of nonlinear 
and/or quasilinear effects, we have made similar com- 
parisons after reducing the impulse amplitudes of the 
enhanced cross-tail potential drop A V ( t ) m our model 
storms 

Figures 7a and 7b show comparisons of Di™ and 
D r ii when the enhanced cross-tail potential drop AV 
for case (a) is reduced by a factors of 2 and 4, respec- 
tively When the average stormtime cross-ti 
tial drop decreases, the agreement between 
Dll' improves considerably The geometric-] 
crepancy amounts to a factor of 1 4 or 1 1, respectively, 
when AV (t) is reduced by a factor of 2 or 4 The agree- 
ment is quite good despite the fact that the quasilin- 
ear diffusion coefficient does not vary smoothly with p 
As the impulse amplitudes in AV(t) are reduced, cor- 
rections to Di L due to resonance-broadening become 


ul poten- 
Dii and 
mean dis- 
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Fig 7 Plots of diffusion coefficients obtained via (14) from tune-reversed (open circles) and time-forward 

(filled circles) simulations, for comparison with the quasilinear diffusion coefficient D (dashed curve) given by 
(16) as an implicit function of p for storms that have the same temporal structure as model storm (a) m Figure 2, 
but with half (left panel) and a quarter (right panel) of the impulse amplitude in A V(t), respectively The average 
enhancements (AV(t)) in cross-tail potential drop over the 3-hr main phase are 90 kV and 45 kV, respectively 
The diffusion coefficients (corrected for resonance-broadening) are represented by the solid curves 
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smaller, as the convergence of the solid curve (D r £ L ) 
toward the dashed curve (Z)^) m Figure 76 shows 
Moreover, the agreement between diffusion coefficients 
obtained from simulations run forward (filled circles) 
and backward (open circles) in time improves The ge- 
ometric mean discrepancy between “time-forward” and 
“time-reversed” diffusion coefficients is only a factor of 
1 2 when AV (t) is reduced by a factor of 2 and only 
a factor of 1 1 when A V(t) is reduced by a factor of 
4 This improved agreement is not surprising since 
time reversal (t — * —t) leaves the quasilinear theory 
diffusion coefficient invariant under time reversal As 
max[AF(t)] — > 0 quasilinear theory increasingly better 
approximates the simulated stormtime transport, and 
so the values of estimated from simulations run 
forward and backward in time become less distinct from 
each other We thus conclude that the remaining dis- 
crepancies between D r A and Df™ are attributable to 
unspecified nonlinear effects of wnich we have not taken 
account 

In a numerical “experiment” we found that excellent 
agreement between £>£*£ and could be achieved 

by arbitrarily increasing Aui/l-n by a factor ~ 3 from 
the value specified by (17). However, we can think of no 
physical rationale for actually postulating such a mag- 
nification of the bandwidth over which DT r should be 
averaged We have considered the possibility that the 
temporal variation of AV(t) in (3) might increase the 
spread m drift frequencies during transport (i e , for 
0 < t < T) beyond that implicit in (17), m which the 
factor (ddz/dL)^ pertains to quiescent drift frequen- 
cies The guiding-center simulations presented earlier, 
however, show little evidence for such an effect Typi- 
cally, the mean drift penod of a representative ion dur- 
ing the mam phase only slightly exceeds the quiescent 
drift period, as it would if the ion were drifting under 
the influence of a constant enhancement AV = (AV (t)) 
of the cross-tail potential drop Thus, the Q.i/2-n which 
appears m (15)-(18) should perhaps have been inter- 
preted as the mean main-phase drift frequency rather 
than as the quiescent one This correction would ap- 
pear to be small, corresponding to a rightward shift 
of the plotted curves in Figure 6 by A/x < 1 MeV/G 
However, the stormtime presence of (AV (<)) does trans- 
form the quiescent drift shell of interest into a storm- 
time band of drift shells over which E(£l 3 / 2 ir) should 
presumably be averaged We will explore the ramifica- 
tions of this refinement elsewhere In our simulations 
the transport-induced spread in quiescent drift frequen- 
cies has in some cases exceeded (17) by ~ 40%, but this 
magnification of Au>/2n would be too little to eliminate 
the remaimng discrepancies ~ 1 9 between D r £ L and 
Dff™ in Figure 6 The very small change m the drift 
frequency during the storm also seems to eliminate the 
possibility of trapped-particle effects A preliminary 
test suggests that the replacement of L 6 m (15) by its 


transport-averaged value (L 6 ) would also be a relatively 
unimportant nonlinear correction 

7 Diffusion Averaged over Storm Ensemble 

It could be argued that quasilinear theory is more 
appropriately applied to an ensemble of model storms 
than to an individual storm We have tested this hy- 
pothesis by randomly choosing 20 different storms hav- 
ing a 184-kV mean cross-tail potential drop from the 
100 storms that we constructed We averaged the dif- 
fusion coefficients obtained from the simulations, stan- 
dard quasilinear theory, and the resonance-broadened 
quasilinear theory over the 20 storms The results are 
shown in Figure 8 [Chen et al , 19926] The^ensemble- 
averaged quasilinear diffusion coefficient D q ^ L (dashed 
curve) and its resonance-broadened counterpart D r A 
(solid curve) are considerably smoother than D^ L ana 
respectively, for an individual storm The ensem- 
ble-averaged diffusion coefficients -Df)™ from the time- 
reversed and time-forward simulations (open circles and 
filled circles, respectively) typically agree much better 
with the theoretical diffusion coefficients m Figure 8 
than does with the theoretical diffusion coeffi- 

cients for an individual storm m Figure 6 The mean 
discrepancy between £>£*£* and is only a factor of 
12 



Fig 8 Ensemble-averaged diffusion coefficients D"™ obtained 
from the time-reversed (open circles) and time-forward (filled cir- 
cles) simulations, for compareon with quasilinear theory (Dy^ 
as a function of p, dashed curve) for equivalent ensembles of 20 
storms, but with impulse amplitudes halved and doubled relative 
to the middle family of curves and data points so as to produce 
the mean A V(t) values shown The solid curve represents the 
ensemble-averaged diffusion coefficients D T £ L corrected for reso- 
nance broadening effects 
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We have also halved the impulse amplitudes of A V in 
the 20 storms of our ensemble and averaged the result- 
ing diffusion coefficients obtained over the 20 storms 
The results are shown via the middle family of curves 
and data points in Figure 8 As expected agreement be- 
tween Z>£*™ and D q LL is even better (geometric-mean 
discrepancy is a factor of 1 1 for {AV(t)) = 92 kV ver- 
sus 1 2 for (AV (t)) = 184 kV ) Corrections of D^ L due 
to resonance broadening are smaller in this case 
When we reduce all the impulse amplitudes in AV (t) 
by a factor of 4, we find remarkably good agreement 
among the ensemble-averaged diffusion coefficients ob- 
tained by the various methods For example, the geo- 
metric-mean discrepancy between D£ l £* and D Y L 
amounts to a factor of 1 03 and the geometric-mean 
discrepancy between the diffusion coefficients deduced 
from time-forward and time-reversed simulations 
amounts to a factor of only 1 02 (see lower family of 
curves and data points, Figure 8) 

8 Summary and Conclusions 


We have used a dynamical guiding-center model to 
investigate energetic charged-particle transport m re- 
sponse to storm-associated impulses in a model of the 
convection electric field Our simple magnetosphenc 
model allows us to compare our numerical results with 
analytical descriptions of particle transport such as the 
quasihnear theory of radial diffusion [Falthammar, 
1965] Thus, we have tested whether quasihnear theory 
can appropriately be applied to describe the charged- 
particle transport caused by electrostatic electric field 
fluctuations over time intervals as short as an individual 
storm Furthermore, we have begun to use our simula- 
tion results to map phase-space distributions from the 
quiet-time to the stormtime nng current by using Li- 
ouville’s theorem A summary of our results follows 
Ions having pi < 3 MeV/G (E £ 110 keV at R = 3) 
gam access to closed drift shells at R ~ 3 mainly by 
direct convection from open (plasmasheet) drift paths 
At L ~ 3 these ions have drift periods that exceed the 
duration of the main phase of the storm The mode 
of access of ions with /i ~ 5-13 MeV/G (E ~ 55-145 
keV at R = 3 appears to be transitional between con- 
vective and diffusive access At L ~ 3 these ions have 
drift periods that are comparable to the length of the 
mam phase of the storm The stormtime transport of 
ions having > 13 MeV/G (E > 145 keV at R = 3) 
resembles radial diffusion across closed drift shells At 


L ~ 3 these ions have drift periods that are smaller 
than the duration of the main phase of the storm. 

The electric spectral density derived from our model 


storm is not very sipooth, and so the quasihnear dif- 
fusion coefficient D?r does not vary smoothly with n 
When we comp area tbe diffusion coefficients Df}™ ob- 
tained from the simulated trajectories with the quasi- 
hnear diffusion coefficient D q £ L {p, L), we nevertheless 


found surprisingly good agreement for 4 distinct indi- 
vidual model storms The aggregate geometric-mean of 
discrepancies between and D q £^ for the 4 storms 
amounted to a factor of 2 2 even though DVr varied ir- 
regularly with \jl by 4-5 orders of magnitude for each of 
the storms studied When we invoked nonlinear drift- 
resonance broadening effects, we found that the dis- 
crepancies between quasihnear theory and D?’™ were 
slightly reduced through a smoothing of the sharp rel- 
ative minima and maxima m The aggregate ge- 

ometric meaij of the remaining discrepancies between 
and for the 4 different storms amounted to a 
factor of 2 0 When we reduced the impulse amplitudes 
in the enhanced cross-tail potential drop A V{t) of our 
model storm, the agreement between D q / L and Df;™ 
improved considerably For example, the geometric- 
mean discrepancy is a factor of 2 3, 1 4, and 1 1, re- 
spectively for a particular storm with (AV) = 180 kV, 
90 kV, and 45 kV This convergence towards 1 0 sug- 
gests that the discrepancies between D s ™ and D^ L are 
attributable to nonlinear effects When we averaged 
Df L and D£™ obtained over an ensemble of 20 ran- 
dom storms, we found even better agreement between 
D 9 j* L and Dfj* (mean discrepancy amounted to a fac- 
tor of 1 2 for (AV (t)) « 184 kV) When we reduced the 
impulse amplitudes in A V(t) for all the storms in the 
ensemble by a factor of 4, we found remarkably good 
agreement between these ensemble-averaged diffusion 
coefficients (mean discrepancy factor was 1 02) 

We developed a simple model of the pre-storm 
(steady-state) phase-space density, obtained by balanc- 
ing radial diffusion against charge exchange, which pro- 
duces features qualitatively similar to those found in 
observations We used our time-reversed simulations to 
map phase space densities / from pre-storm to storm- 
time in accordance with Liouville’s theorem The 
stormtime transport of ions having /x £ 13 MeV/G (E 
> 145 keV at R = 3), for which the transport resembles 
radial diffusion, seems to produce little change m dnft- 
averaged phase space density /. However, the storm- 
time transport of / produced a major enhancement in 
/ from the pre-storm phase-space density at energies 
~ 30 — 145 keV, which are representative of the storm- 
time ring current These are energies for which our sim- 
ulations have shown that many of the ions are trans- 
ported on the mghtside from open trajectories to the 
final drift shell (L ~ 3) of interest 

Our preliminary results on the mapping of phase 
space densities are particularly satisfying since they re- 
produce many observed features of ring-current phase- 
space distributions Thus, we are performing addi- 
tional mappings to other Ir-shells of interest We plan 
to refine our treatment of the pre-storm and boundary 
phase-space distributions so as to achieve a more real- 
istic model We also will include loss processes such as 
charge exchange in our simulations so that we can map 
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phase space distributions through the recovery phase as 
well as through the main phase of a model geomagnetic 
storm 
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